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I. INTRODUCTION 



The Kohn-Sham (KS) formulation 1 of Density Functional Theory (DFT)2~- is one of the 
most important tools for the calculation of electronic structure of molecules and solids. In all 
practical applications of DFT, however, approximations to the exact functionals have to be 
made* — . Exact relations for density functionals and density functional derivatives can play 
an important role in the development of accurate approximations to the exact functional*^. 
In this paper the following equation satisfied by the exact correlation energy and potential 
is derived: 



Here p^ and p^_ x are the ground state charge densities of an interacting N and (N — 1) 
electron system of the same Hamiltonian with multiplicative external potential i>g Xt ([pat]). 
The potential i>g xt ([pjv]) is constructed to keep the charge density of the N electron system 
independent of the coupling strength parameter y^~— that scales the electron-electron in- 
teraction strength. At 7 = 1 full strength Coulomb interaction between electrons is included 
and the external potential i>g Xt ([pjv]) is the external potential of the fully interacting system, 
while 7 = corresponds to the non-interacting Kohn-Sham potential. 

It is well known that the correlation energy EJ [pn] at coupling strength 7 can be ex- 
pressed as the integral^— 



J 

where V^ e [pjv] is the electron-electron Coulomb interaction energy and Eh x [pn] is the sum 
of the Hartree and exchange energy. An alternative way of expressing this relationship is 



This equation shows that the there is a close relationship between the correlation energy 
El \pn\ at 7 and the correlation energy and the correlation part of the kinetic energy T 7 \pn\ 
for all coupling strengths between and 7. This expression is not solely in terms of the 
correlation energy. Equation (JT]), on the other hand, provides a relation entirely in terms of 



E" c [ P n] - El J 




(1) 




(2) 
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correlation functionals. It is remarkable in that it not only relates the correlation functionals 
at different densities of the same Hamiltonian, but indicates a relationship between the 
correlation functionals at all coupling strengths between and 7. Equation ([T]) can serve as 
a test for approximate correlation functionals provided accurate densities are used. 

II. PROOF 

In the adiabatic connection approach^"— of the constrained minimization formulation of 
density functional theoryi^^^ the Hamiltonian H 1 for a system of N electrons is given by 

W = f + 1 V e e + vl t [ PN }. (4) 

Atomic units, h = e = m = 1 are used throughout. T is the kinetic energy operator, 

1 N 

r = 4£ v * 2 ' (5) 

i=l 

and 'yVee is the scaled electron-electron interaction operator 

N 1 

• ^ ■ \ i i\ 

The external potential 

N 

Cxt M = v ^t ([pjv] ; r*) , (7) 
i=i 

is constructed to keep the charge density of the iV-electron system fixed at pn (r) , the ground 
state charge density of the fully interacting system (7 = 1), for all values of the coupling 
constant 7. The external potential has the forrn^ 1 ^ 

*4t(M ; r) = (1 - 7) v hx ([p N }; r) 

+ vK\Pn]; r) - v2([p N ]; r)+vl t ([p N ] ; r), (8) 

where vl xt ([pN-} ; r) = v cxt (r) is the external potential at full coupling strength, 7 = 1, 
and Uextd/'Ar] > r ) ^ s the non-interacting Kohn-Sham potential. The exchange plus Hartree 
potential^ Vh x {[pN}', r), is independent of 7, while the correlation potential v2([pn]',*) de- 
pends in the scaling parameter 7. 
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The energy functional F 7 [p] is defined as^ 1 ^ the sum of the kinetic and interaction 
energy, 



F 7 [p] 



v[r7 
p 



min ( \I> 



T + -yV e 
f + >yV € 



p 



(9) 



where according to the Levy constrained minimization definition^, the wavefunction 



yields the density p and minimizes ( \l/ 



$ ) . For unrepresentable densities^ 1 ^, the 



T + 7I4 

densities are derived from the groundstate eigenfunctions of a Hamiltonian: 

^|*?r> = ^|*M 

^ = f + 7 y ee + c t W • 



(10) 



Note that by construction of i>e Xt {[pn] ] r ) > Eq. flHl) Pat = Ptv> is independent of 7, but the 
groundstate density of the (N — l)-electron system p^_ 1; is a function of 7. -F 7 [p] can be 
decomposed as^ 

F>] = T>] + 7 £ te [p] + £J [p] , (11) 



The correlation energy EJ [p] is defined asM 



E7 [p] 



^p 



p 



-(* 



T + -yV e( 



(12) 



where \^° p ) is the Kohn-Sham independent particle groundstate wavefunction that yields 
the same density as the interacting system at coupling strength 7. Eh x [p] is the sum of the 
Hartree and exchange energy 



E hx [p] = (# 



Ve, 



(13) 



and the non-interacting kinetic energy functional T° [p] is given by 



T°[p] = (K 



T 



(14) 



The full kinetic energy 



t 1 [ P ] = h 



T 



P 



T° [p] + V [p] , 



(15) 



with the correlation part of the kinetic 



T 



T 



(16) 



Assuming that i^ 7 [p] is defined for non- integer electrons^ 1 ^, at the solution point 

5F~< [p] 



5p(r) 



+ <xt(b]; r ) = p 



(17) 



where p is the chemical potential. 

In a recent paper— the author proved that the correlation part of the kinetic energy 
satisfies 



V [pn] - T2 [pl_ x ] 

d *r>(p N (0-^(0)™ 



(18) 



In appendix A, Eq. (lAlOp . it is shown that 



T2 [pI^] 

= m [pi-*] - ^ \pi_ x \ 

+7 J l ^r dpN ~ 1 (v2 ([p N ] ; r) + -fv hx {[p N ] ; r) - jv hx ([p^_J ; r)) , (19) 
where {[p] ; r) = j^Fj E c [p] and v hx {[pn] \ r) = j^E hx [p). It follows that 

d 2 

+ j ^ ^Pn^I ( r ) ^2 Qp^] - r ) + 7 ^ Q pjv ] - r ) _ 1Vhx ( ^7 _J . r ^ 

+7^ y dV ^ 1 ' (v2 ([p N ] ; r) + jv hx ([p N ] ; r) - 1Vhx ([pi _J ; r)) . (20) 

A similar expression can be derived for -§^T^ [pjy] but since pat is independent of 7, the result 
is simply 

fa m = -j-^m k_j . (21) 

Taking the derivative of (fT9|) with respect to 7, using ( 120]) and ( )2T|) and some algebra, it can 
be shown that 



d 2 „„, r , d 2 



dl 2 c Lr.vj 

+ J d3r ^PN-i ( r ) ( v ux ([Pn] ; r) - v ux ( ; r)) 

£ / rfV ^< r '»- *-> < r '»G^W^ w ) (22 > 



7=0 



This equation can be integrated with respect to 7 and since by definition E® [p] = J^E] [p] 
0, 

fa M - fa 

r <9p 7 (r) 

+ y rf3r ^^ ( v c (M ' r ) + TVux ([Pn] ; r) - 7^ ( [p^_J ; r)) 

From Eq. (1A10I) applied the 7 invariant iV-electron density, it follows that 

^ M = MM , (24) 



Substitute f ETOjl and (J23D in with the result 



£ c 7 [Pat] - £ c 7 [p^_J 

+7 y d 7 J d 3 r— p 7 N _ t (r) ([p N ] ; r) - u ux ( [p^J ; r)) 

= I dV ( Pjv (r') - pl_ x (r')) M • (25) 

This equation still contains reference to the Hartree and exchange potentials. In order to 
find an equation that explicitly only contains correlation energy functionals and potentials 
consider the following. In a recent paper— the author also showed that the full kinetic energy 
functional satisfies 

[p N ] ~ T 1 [pl_ x ] 

= /^(p,M-^-i(rO)^ (26) 



6 



Taking into account that 

T 1 [p] = T° [p] + 77 [p] (27) 
and that pjj_i is a function of 7^, the 7 derivative of f|26|) can be written as 



3 ( r dV ( PN (r') - pl_ x (r')) M 



, dp^ (r) / ST° [p^] 5T q [p N ] 



+ d^r-^^—— / f J - . (28) 



Comparison with ffTBj) shows that 



By construction, from (TTT|) and (TTTT) 
<5T° [ PJV ] 



<W (r) 

ST° [p N ] 



+ vl ([p N ] ; r) + jv ux ([p N ] ; r) + v^ t ([p N ] ] r) = //^ 



v , • v , + ^ 7 ( [Pn-i] ; r ) + TVux ( Hr-J ; r) + w e 7 xt ( [p^_J ; r) = p? N _ x (30) 
and hence 

^Pat-i ( r ) <W (r) 

= ( K-J I r ) + l v ux ( [Pn-i] ;r) -v2 {[p N ] ; r) - ~/v ux {[p N ] ;r) + p? N - ff N _ v (31) 



Now 



£("-l)=/*r*%*U. (32, 



From (I32l) . (l3Tj) . ( ]29l and (|25|) the main result of this paper follows 



E2 [p N ] - E2 [p 7 _J 



d 3 r'(p N (v')-pl_ 1 (v'))v c ([p N };r). (33) 
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III. DISCUSSION AND SUMMARY 

The identities (lisp and were derived for unrepresentable densities^ and the proof 

of Eq. flAlOp also makes use of the eigenfunctions of a Hamiltonian. Therefore Eq. ([T]) is 
valid for densities derived from groundstate wavefunctions of a many-particle Hamiltonian. 
The assumption was made that all functional derivatives are well behaved and this implies 
that the functionals are defined for non- integer particle numbers^ 8 . Equation ( ([T|) is in- 
teresting in that it provides an expression entirely in terms of correlation functionals and 
functional derivatives. It is valid for exact densities and any test of an approximate correla- 
tion functional will only be reliable if exact densities are used. The 7-dependent correlation 
energy can be derived for any approximate correlation functional since^ 



where p\ (r) = A 3 p (Ar) is the uniformly scaled density. However, to implement ([T]) accurate 
densities for the (N — l)-electron system is needed for all electron interaction strengths 
between zero and 7, and this may be difficult to implement in practice. 

In summary, a new identity for the correlation functional was derived that relates the 
functional and functional derivatives of different densities of the same Hamiltonian but where 
the particle number differs by one. 



Appendix A 



El [p] = ^K' Pi 
. 7 



(34) 



FromEq. (JT2D 



E2 [p£_J - T c 



1 




(Al) 
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The derivative of EJ [pJv-i] with respect to 7, from definition ffT2]) . can therefore be expressed 



as 



7 



+ <* 7 



T 



P + 7^ 



^ 7 



"7 



— * 7 7 

Qsy Pjv-1 



Vl> 

Pn-i 



d_ 

dj ' Pn- 



(A2) 



Upon adding and subtracting (c.c. stands for the complex conjugate of the previous term) 

d 



+ 



— # 7 7 

C?7 Pjv-i 
* 7 

^7 p iv-i 



'iV- 



-l,extM * 



+ C.C 



7=0 



l^l.cxtMl^p^ )+CC 



and utilizing the normalization of the wavefunctions which implies that 







Q'y \ Pn-i Pn 



(A3) 



(A4) 



Eq. (ITbl) becomes 
^7 



£ c 7 K-i] 



£ c 7 K-i] - r c 7 K-i] 
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Pn- 



d 



+ < -^pX., ICt [pat] I *° ol , > + ( *° ol , ICt M| ^ 



dj 



Pn-i 



Qru P N „ 

d 



Pn-i 



Pn-i 



v;, 



\ ~ Pk-1 

This equation can be simplified since 
d 



vl> 

7 

Pn-i 



(A5) 



— ^ 7 7 
c?7 p jv-i 



^ 7 



ICxt [p 



TV 



\1> 7 7 

Qrv Pn-1 



d 3 r 



dpl^ (r) 
97 



*c 7 xt ([Piv] ;r) 



(A6) 



From definition of the exchange energy^ 



Pn- 



v P , 



Pn-i 



E x [pl_ x ] + U [pl_ x ] 
Ehx [Pn-i] , 



(A7) 



is the sum of the exchange E x [p^-_ 1 ] and Hartree interaction energy U [pjv_i] of the (N — 1) 
electron system. The charge density p y N _ 1 is a function of 7^, therefore 



Qry \ P N -1 



d 6 r 



dPN-i ( r ) 
(?7 



Wfe {[Pn-i] ;r) , 



where 



( [pX--J ; r ) 



(A8) 



(A9) 



^-1 ( r ) 

is the sum of the exchange and Hartree potentials for the (N — l)-electron system. Using 

yield the same density p]^_i, 



Eqs. (TMD, 



and the fact that 



\l> 7 7 ) and 

Pn-i 



Pn- 



Eq. (IA50 can be expressed as 




97 



£ c 7 K-i] - T c 7 K-i] 



7 



$7 



(M I r ) + 7^te (M ; r) - ~/v hx ([p^_J ; r)) (A10) 
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